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l Introduction 

This report covers wort carried out on NASA contract NAS8-36955 during the year 1990. The report is divided 
into the following sections :- 

Preionisation Review:- A literature survey to identify the reguired parameters for effective preioni- 
sation in TEA CO* lasers and the methods and techniques for characterising preionisers. 

Cavity Rodel:- A numerical model of the LP-140 cavity used to determine the cause of the transverse 
node stability improvement obtained when the cavity was lengthened. 

LP-140 Discharge Characteristics:- The measurement of the voltage and current discharge pulses on the 
LP-140 and their subsequent analysis resulting in an explanation for the low efficiency of the laser. An 
assortment of items relating to the development of high-voltage power supplies is also provided. These items 
are guidelines that were prepared for a GSRP student. 

Frequency Chirp:- A program for analysing the frequency chirp data files obtained with the HP time 
and frequency analyser. A program to calculate the theoretical LIHP chirp is also included and a comparison 
between experiment and theory is made. 

Laser Linewidths:- A program for calculating the CO* linewidth and its dependence on gas composition 
and pressure. The program also calculates the number of axial modes under the FWHH of the line for a given 
resonator length. A graphical plot of the results is provided. 

Appendix (1):- This is a listing of the basic program used to extract the energy deposition rate, 
inpeadance etc. from the discharge voltage and current pulses. 

Appendix (2):- Two Nathcad listings, one to calculate arbitrary TEHi,» nodes up to 3,3 and the other 
to calculate the focal length of a thick lens given the radius of curvature of the surfaces. The second 
document is useful for converting ROC to focal length and vice-versa for the resonator code in section 3. 

Appendix (3):- This is a copy of a joint paper between Dr. Jaenisch and myself to be presented at the 
SPIE High Energy Lasers Conference in Los Angeles on January 24 th 1991. 

Due to the variety of topics covered an overall conclusion is not provided as each section is self- 
contained. 
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2 Preionisation Review 

2 . 1 Introduction 

This review is a copy of a document submitted previously, only minor formatting changes have been 
undertaken. The object of the review was to provide a brief overview of the means of obtaining quantita- 
tive data on the relative merits of various preionisation schemes. 

The discharge in a TEA CO* laser can be divided into three phases, firstly preionisation, which involves 
the injection of electrons into the gap in which the discharge will be formed. Secondly the application 
of a high dV/dt voltage pulse to avalanche the preionisation electrons up to an electron concentration 
suitable for maintaining a stable discharge. The third phase consists of the deposition of energy into 
the discharge. During this final stage the voltage across the electrodes collapses to a value essen- 
tially independent of the driving circuit. Allthough the preionisation pulse can be completed prior to 
the application of the avalanching pulse, it is more common to arrange the timing such that the 
avalanching occurs during the peak of the preionisation pulse. This ensures that a minimum number of 
preionisation electrons are lost to attachment and also limits the degree of avalanching required to 
form the discharge. The preionisation pulse usually does not continue during the third phase, except for 
long pulse operation where electrons must be continuously injected into the discharge to maintain its 
stability. For large cross-section discharges (-10x10 cm 2 ) the preionisation can not be generated uni- 
formly throughout the discharge gap and so to overcome this, a high preionisation density is created 
adjacent to one of the electrodes and a weak (non-avalanching) voltage applied across the electrodes to 
drift this preionisation uniformly through the discharge volume prior to the application of the ava- 
lanche voltage. The preionisation of the discharge volume is critical to the successful operation of the 
discharge and is therefore a vital consideration in the design of a transversely excited gas laser. 

2.2 Preionisation requirements 

The main requirements of a preionisation source are:- 


2.2.1 Large electron production 



quantitative agreement of these values. The actual discharge has an electron density of 10 12 -10 14 
cn 3 and experimental evidence' 2 5 ’ has shown that initial preionisation densities 

above the minimum required value will result in more efficient energy extraction from the discharge 
volume. Additionally with a larger initial electron density less avalanching is required to obtain a 
stable discharge and a lower initial voltage can be applied to the electrodes. 

2.2.2 Uniform electron production 

The preionisation must be produced uniformly throughout the discharge to prevent regions of high 
current density from occur ing. This is important as these high current density regions can readily 
degenerate into arcs. 

2.2.3 Efficient electron production 

Within the requirements of the LAWS project, only a limited power budget is available and so the 
efficiency of the preionisation process is important, therefore the number of electrons generated 
per unit energy input to the preioniser should be as hi#i as possible. 

2.2.4 Contaminant free preionisation 

The lifetime requirements of the LAWS program requires that the preionisation source should provide 
little or no contamination that nay limit the lifetime of the laser. As an example spark preionisa- 
tion has in the past been preferred for flowing-gas lasers as it provides large numbers of elec- 
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trons. However in a sealed laser' 2 ' 9 ' 2 ' 10 ’ the high temperature spark results in alnost complete 
dissociation of the COa in the region of the spark resulting in a shorter lifetime for the laser 
than if a diffuse discharge preionisation source had been used. Additionally the spark pins have a 
tendency to sputter material onto the main discharge electrodes and the laser optics. 


2.3 Preioniser test methods 

Traditionally most preionisers have been tested by placing them in a laser and provided the laser has 
worked the preioniser has been considered satisfactory. There have been very few attempts (in the open 
literature) to compare preionisers directly. The techniques that have been used to examine the charac- 
teristics of preionisers are:- 


2.3.1 UHV chamber 

A OHV chamber' 2 ' 6 ' 2 ' 9 ' 215 ’ chamber enables the direct measurement of electron production from the 
preioniser to be undertaken. Typically a preioniser is placed inside a chamber such that the 
electrons generated are collected by a very weak electric field across a pair of electrodes. These 
collection electrodes can be moved to measure the preionisation from various portions of the preio- 
niser. The dependence of the preionisation density on the preioniser drive circuit parameters can 
obviously be determined using this technique. The signal obtained by this technique for low electron 
densities is very small (*nV) and considerable care is required in the design of the apparatus. A 
dummy electrode is utilised to allow a differential signal to be obtained. Additionally suitable 
filters are normally used to eliminate HP noise which arises from the preionisation discharge. As it 
has been shown ' ' ' that small quantities of organic contaminants can greatly 

increase the preionisation density for uv based preionisers, to eliminate fluctuations in measure- 
ments extreme care must be taken to ensure the cleanliness of the chamber is maintained. By repeti- 
tively pulsing the preioniser in the chamber the dependence of the preionisation density on 
preioniser age could also be determined and after many pulses the gas in the chamber could be 
analysed to determine the degree of degradation due solely to the preioniser under test. However 
with the addition of life testing of the preionisers, a high pulse rate becomes desirable to allow 
the tests to be completed in a reasonable time. This would require cooling and gas circulation 
within the chamber. From the literature the effect of the preioniser on the gas mixture can easily 
be determined after *100,000 pulses (for spark preionisers' 2 -*' 2 - 10 ’). if an hour is allowed for 
this test, this provides for a pulse rate of 28 Hz, which is minimal. As the preionisers are intended 
to be very clean a pulse rate of 50 Hz nay be more suitable. 


2.3.2 Photographic studies 

Both time-averaged' 216 ' 2 ' 20 ’ and time resolved' 2 ' 16 ' 2 ' 21 ' 2 ' 22 ’ photographic studies have been used 
to study preionisation and discharge formation. The time averaged studies although of limited value 
will show regions of very high current density. Time resolved photography provides much more infor- 
mation and can be used to follow the discharge development and enables flaws in the uniformity to be 
easily seen. The technique requires the use of a streak camera or preferably a framing camera with a 
resolution -10-20 nS. 


2.3.3 Light: Emission 

Many of the preionisation schemes rely on the production of uv light which is used to preionise the 
discharge volume. Measurements have been made of the light production from these preionisers using a 
photodiode * , however the disadvantage of this technique is that the correlation between light 

output and electron density, although probably valid, has not been demonstrated. 
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2.3.4 Microwave interferometry 

Early atteipts to measure the electron concentrations in a discharge used an X-band microwave 
interferometer' 2 ' 23 however the sensitivity of the device was insufficient and the technique would 
appear to have been abandoned. 

Of the techniques outlined above the UHV chamber is the lost important for the quantitative coiparison 
of preionisers. The tiie averaged photography is relatively easy to accomplish, although ideally a large 
format camera should be used to ensure adequate resolution. The time resolved photography provides the 
only means of studying the discharge formation on the preioniser thereby enabling transient 'hot spots', 
which frequently indicate potential failure locations, to be detected. The microwave interferometer can 
be discounted due to its limited sensitivity. 

2.4 List of required equipment 

The equipment ideally required for the characterisation of the preionisers therefore consists of the 
following items. It should be noted that only a general description is given at this point. 

2.4.1 UHV chamber 

An UHV chamber assembly consisting of collection electrodes, ports for supplying gas and to enable 
gas sampling. Suitable feedthroughs for the electrical connections will be required. A large port to 
which the test preioniser assembly can be fixed and a large transparent window directly opposite 
this port to enable photography of the preionisation discharge. In order to facilitate movement of 
the collection electrodes, rotary feedthroughs must be provided if manual adjustment is envisaged. 
Alternatively and possibly preferably if automatic data collection is required a pair of stepping 
motors and their associated electronic drivers is required. If life testing is envisaged then gas 
cooling and recirculation fans must also be supplied. 

2.4.2 Drift field power supply 

The drift voltage required obviously depends on the electrode spacing but typically a voltage of 
-0.3 Vci'Horr" 1 is used giving a drift voltage of -2 kV for an 8 cm gap at atmospheric pressure. 
Thus a 0-5 kV dc power supply should suffice. 

2.4.3 Detection circuitry 

This consists primarily of a differential ac amplifier and HF filters (which may be included in some 
amplifiers). Traditionally the output from the amplifier has been fed to an oscilloscope and dis- 
played so that a photograph of the signal can be obtained. This obviously lends itself to digitisa- 
tion. 

2.4.4 Gas analysis 

If the effect of the preionisers on the gas chemistry is to be studied a mass spectrometer will be 
required. 

2.4.5 Time integrated photographic analysis 

A suitable large format camera is required. 

2.4.6 Time resolved photographic analysis 

A suitable high speed framing camera is required. 
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The above list only covers the equipment required to analyse the electron production from the various 
preionisers. Thus in addition to this list suitable circuitry for driving the respective preionisers 
would be required together with high voltage and high current probes for their electrical analysis. 
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3 Cavity Nodal 

3.1 Introduction 

The transverse node structure of the laser output pulse is important for obtaining a good far-field beai 
profile. The initial bean profile of the LP-140 laser had considerable transverse node instability <il) . 
A schematic of the LP-140 cavity is shown in figure 3.1. Jaenisch <31) found that when the cavity length 
was extended by increasing the seperation between the grating and the lens, the transverse node stabil- 
ity inproved. The nodel discussed here was developed to provide an explanation for this inproveient. 



Figure 3.1. The LP-140 cavity layout. 


3.2 The matrix model 

The nodel is based on the work of Siegnann (32> ( and has been inplenented as a fathead* 33 ’ docunent. In 
this nodel each of the optical elenents in the cavity can be expressed as a 2x2 natrix. By conbining the 
natrices a single natrix for the overall optical systen can be obtained. Fron this natrix the nagnifica- 
tion and Fresnel nunbers of the resonator can be deternined. 

The sinple geonetrical nagnification of the resonator varies with the conplex curvature, g of the input 
wavefront and for lasing action to occur the wavefront mist be capable of reproducing itself over nany 
round trips. This leans that the conplex curvature of the wavefront after one conplete round trip lust 
return to its initial value. This condition constrains the nagnification to a single value and it is 
this value that is calculated by the prograi. 

Rather than provide a detailed theoretical treatise of the m>del, the reader is referred to the 
book* 32 ’ by Siegnann which provides a conprehensive account of the theory. The fathead docunent 
follows :- 
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LP-140 CAVITY MODEL 


This document uses the matrix method to model the LP-140 cavity. Each of the 
cavity elements is represented by a matrix where the variables for each matrix 
are : - 


Output coupler radius of curvature (m) 

Cavity length between output coupler and lens (m> 
Radius of curvature of first surface of lens (m) 
Refractive index of lens material 
Thickness of lens material (m) 

Radius of the second surface of the lens (m) 
Spacing between the lens and grating (m) 
Wavelength of laser radiation (m) 

Mirror diameter (m) 

Unstable resonator output coupler half width (m) 
Total cavity length (m) 


R1 

= -5 

D1 

= 2.40 

R2 

= -14.05 

N1 

= 2.407 

T1 

= 0.00635 

R3 

= 1000000 

D2 

= 0.1 

X 

= 0.0000106 

W 

= 0.04 

a 

= 0.009525 

LEN 

= D1 + D2 + Tl-Nl 


This provides the following component matrices for one beam direction 


Output Coupler Cavity Space 


First surface of lens 



Lens thickness Second lens surface Grating Spacing Grating 



For the beam travelling in the opposite direction we must reverse the sign of 
the curvatures on the lens surfaces to give:- 


Grating side surface Cavity surface 


R4 := -R3 R5 := -R2 



The total system matrix is obtained by multiplying out the matrices :- 

M := M M -M M -M -M -M M -M -M -M M 
s 12345676 52 4 32 2 

Giving: - 



0.504 3.765 

M 

= 1 

(This is a check on 

M 

= [0.003 2.01 

s 


the validity of the 


s system matrix and 

should equal 1 . ) 



The ABCD elements of the system matrix are:- 

A := M B := M C := M D := M 

s s s s 

0,0 0,1 1,0 1,1 

The half trace parameter, m is given by:- 

A + D 

ra := 

2 

giving m = 1.257 (For -1 < m > 1 the resonator is unstable, whilst for 

-1 < m < 1 the resonator is stable.) 

If the resonator is unstable the magnification is given by:- 




M := if 


m > 1 ,m + 


m. 



M 




CM 

if 

m < -1,-m - N 

ro - 1 , M j 


i.e. M = 2 . 018 ( If the magnification has the value 0 this indicates the 

stable resonator case.) 


Now the Fresnel number is calculated for the resonator, together with the 
collimated and equivalent Fresnel numbers if the resonator is unstable. 


Fresnel Number 


2 

W 

N := 

f 4 • LEN • X 


If unstable then:- 
Collimated Fresnel Number 


N := if 
c 


M- a 

M ji 0, ,0 

B-X 


and Equivalent Fresnel Number 


N 

eq 


:= if 


M / 


M 


2 


0 , 



2 

2 • M 


,0 
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The various Fresnel numbers are:- 


N = 15.003 N = 4.587 N = 1.73 

f c eq 

where a value of zero for the collimated and equivalent Fresnel numbers 
indicates a stable resonator. These parameters together with the magnification 
enable a preliminary estimation of the resonator mode stability and losses to 
be obtained. 
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3.3 Analysis 

The oodel was run using the initial cavity configuration and then for the cavity as nodified by Jae- 
nisch <31) . The iodel was then rerun with the initial cavity length but with a nodified prescription for 
the lens, LI (figure 3.1). The object of this exercise was to see if the shorter cavity length could be 
maintained whilst also retaining the magnification and equivalent fresnel numbers of the extended cavity 
configuration. It should be noted that the output coupler, HI (figure 3.1) could also have oeen nodified 
to attempt this, however due to the high cost of modifying this element, this option was not undertaken. 

A summary of the results from each of the models is presented in table 3.1. It was found that when the 
prescription of LI was varied only one of the desired properties could be matched for a given prescrip- 
tion. Thus two models appear in the table, one to match the equivalent fresnel number and the other to 
match the cavity magnification. 


Cavity 

Cavity 

Lens, L2 

Equivalent Fresnel 

Magnification 

Type 

Length (n) 

ROC (■) 

Number, N. q 

H 

Original 

2.50 

14.05 

1.73 

2.02 

Extended 

2.96 

14.05 

1.63 

2.11 

Model 1 

2.50 

14.80 

1.84 

2.12 

Hodel 2 

2.50 

13.48 

1.63 

1.93 


Table 3.1. Cavity Modelling Results. 


The existing lens is a plano-convex ZnSe element with the convex side facing the discharge. Only the 
radius of curvature (ROC) of the convex surface was varied. It can be seen from the table that a radius 
of curvature of 14.8 m is required to match the magnification whilst a radius of curvature of 13.48 m is 
required to match the equivalent Fresnel number. Siegmann <3a ’ has shown that the optimum value of the 
equivalent Fresnel number for maximum transverse mode stability is (n + 0.5) for a circular mirror 
resonator and (n + 0.375) for a strip resonator, where n is an integer. In our case with a rectangular 
beam cross-section, the strip resonator condition applies. It is obvious that increasing the cavity 
length improves the transverse mode structure by some means other than changing the equivalent Fresnel 
number as this remains approximately constant for both the original and extended cavities. It is possi- 
ble that the unstable transverse mode structure is due to parasitic oscillations from the electrodes. 
Increasing the cavity length results in an increase of the cavity diffraction losses and the elimination 
of the parasitic oscillations. When the ROC of the lens is changed to match the extended cavity magnifi- 
cation, the diffraction loss required can be obtained without the need to extend the cavity. 

3.4 References 

3.1) Jaenisch, H.K.; "Polarisation effects Tasks 3 and 4: LP140"; Final Report NASA Contract No. NAS8- 
36955; Delivery Order 33; (August 1990) 

3.2) Siegmann, A.E.; "Lasers"; University Science Books; Hill Valley, CA.; ch. 14-23 (1986) 

3.3) Hathcad Software Program Version 2.5; Hathsoft Inc.; Cambridge, HA 
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4 LP— 140 Discharge Characteristics 

4.1 Introduction 

The discharge circuit is an inportant component of the overall laser design as it can dramatically 
effect the efficiency of the laser. The circuit is designed to perfori several functions. 

4.1.1 Discharge Initiation 

In order for a stable discharge to be fomed the initial gas voluie between the electrodes mist 
undergo electrical breakdown. For this to occur boiogeneously, preionisation electrons mist be 
injected into the discharge gap prior to the sain voltage pulse appearing across the electrodes. 
Hormally a preionisation source driven by the discharge circuit provides these initial electrons. 
Following preionisation, the sain voltage pulse is applied across the lain electrodes to produce 
breakdown. Breakdown occurs through the rapid avalanching of the preionisation electrons to produce 
a plasia. For this process to occur successfully the lain voltage pulse lust have a fast risetiie. 
The tiling between the peak of the preionisation electron pulse and the initiation of the avalanch- 
ing voltage pulse mist be optiiised to utilise as uany of the preionisation electrons as possible so 
that a minimum maber are lost to electron attachient processes. Much of this tiling can be arranged 
to occur automatically by a suitable design of the discharge circuit. 

4.1.2 Energy deposition 

once the discharge has been formed, it should be optiiised to efficiently pump the upper vibrational 
level of the C0 2 lasing transition and other associated levels (such as the Na vibrational levels). 
It has been shown that there is an optiiui electric field strength per nolecule, known as the 
reduced electric field E/N for pumping the levels of interest. Correct matching of the discharge 
circuit impeadance to the discharge iipeadance, together with the correct choice of voltage for the 
discharge circuit pulse foning network (pfn) ensures this criteria is met. 

4.1.3 Pulse termination 

If the circuit has been correctly designed to the above criteria, the discharge will normally 
extinguish itself homogeneously when the energy stored in the pfn has been deposited into the 
discharge. Occasionally for pulse forming reasons, a tail-biter switching technique may be used to 
terminate the discharge prior to the complete discharge of the pfn. 

4.1.4 LP-140 Discharge circuit 

The LP-140 discharge circuit is shown in figure 4.1. The circuit is duplicated for each of the eight 
discharge electrode pairs, with the exception of the thyratron and 880 HD bypass resistor which are 
common to all the circuits. There are essentially two discharge circuits in this circuit, the first 
is the main discharge circuit consisting of the 0.1 nF capacitor and the second comprises the three 
22 nF capacitors. As there is an inductor in series with the 0.1 nF capacitor, this circuit will 
have a large time constant with respect to the 22 nF discharge circuit. When the thyratron is 
triggered, the 22 nF circuit provides a low energy fast risetime voltage pulse across the laser 
head. The pulse breaks down the gap between the grid and the pin electrodes to provide a preliminary 
preionising discharge. When the slower main discharge pulse arrives at the laser head, the main 
discharge gap is preionised and homogeneous breakdown occurs. The small 500 pF capacitor prevents 
energy from the slower main discharge pulse from being deposited into the preionisation discharge. 
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Figure 4.1. The LP-140 Discharge Circuit. 

4.2 Experimental measurements 

The voltage and current characteristics of the LP-140 were aeasured to enable the discharge to be 
characterised. 

4.2.1 Voltage and current pulse measurements 

The voltage pulse was aeasured using a Tektronix P6015 high voltage probe connected to a LeCroy 9450 
dual 350 KHz digitising oscilloscope. The discharge current pulse was leasured sinultaneously using 
a Pearson Electronics Hodel 110 induction coil. After being digitised by the oscilloscope, the 
pulses were either ported to a plotter or to an IBM coipatible conputer for analysis. Figures 4.2 
and 4.3 show representative voltage and current pulses froi one of the electrode pairs. 

4.2.2 Discharge pulse analysis 

It can be seen fron the voltage pulse that there is an initial sharp voltage pulse rising to 
-4.5 kV. This is the preionisation pulse and it is ianediately followed by the nain discharge pulse 
which rises to -5.5 kV. The voltage pulse then proceeds to ring across the discharge with an 
approxinate average value of -4.5 kV. This ringing, allthough undesirable is not excessive. The 
current pulse rise is delayed with respect to the voltage rise. This is to be expected as the 
discharge inpeadance is required to fall before a significant current flow can be obtained. The 
current peaks at a value of -650 A as the voltage pulse starts to fall. 

Further analysis of the discharge was carried out to enable the discharge iapeadance, electron 
concentration, electron drift velocities and energy deposition to be deterained as a function of 
tine. 

The basic prograa Anlysvi listed in Appendix 1. was used to obtain these characteristics fron the 
initial data files. 
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Figure 4.2. The Discharge Voltage Pulse. 
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Figure 4.3. The Discharge Current Pulse. 
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4. 2. 2.1 Energy Deposition 

The energy deposited into the discharge is given by:- 

E = f l'x/dt (4.1) 

where E is the pulse energy, V and I are the voltage and current respectively and t is the tiie. 
This can be taken as:- 

f= ^(Kx/x6() (4.2) 

n 

where n is the nuiber of digitisation intervals over the discharge pulse and bt is the saipling 
period. Using this enables the prograi to calculate the energy deposition. The energy deposition 
for the voltage and current pulses given earlier is shown in Figure 4.4. 


Energy 

(J) 



Figure 4.4. Energy Deposition in the Discharge. 

It can be seen that a total of -3.4 J is deposited into the discharge. This provides a total of 
27.2 ± 3 J deposited into all eight pairs of electrodes. For a noninal optical output of 1 J 
this provides an electrical to optical conversion efficiency of 3.7 ± 0.2 This value is low 
coipared to a potential efficiency >10 \. Only the energy deposited into the C0 2 (001) and the 
Mv=l-8) vibrational energy levels is useful for lasing action. Lowke et al. <41> have shown 
that the optima excitation of these vibrational levels occurs for an E/N of -(l-3)xl0‘ ls Vci 2 . 
The value of E/N is obtained fron the voltage by:- 


E/N = 


V 

d'x N 


(4.3) 
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where d is the discharge electrode seperation and N is the nolecular nuiber density given by:- 


PxVoL* X A 
RXT 


(4.4) 


where P is the total gas pressure, Vol is the discharge volune, NA is Avogadro's constant, R is 
the nolar gas constant and T is the gas tenperature. The range of E/N values during the 
discharge pulse is divided into intervals of lxlO" 16 Vet' 2 and the energy deposited into each of 
the intervals is calculated (figure 4.5). 



Figure 4.5. The Energy Deposition Dependence on E/H. 

It can be seen that nost of the energy is deposited in the (4-6)xl0 _ls VcT J region and this 
accounts for the low electrical to optical conversion efficiency obtained. To inprove the 
efficiency, the discharge circuit inpeadance should be optiiised to prevent the ringing seen on 
the voltage pulse. Ideally the circuit inpeadance should natch that of the discharge. The 
discharge inpeadance can be approxinated by dividing the voltage by the current (figure 4.6). It 
can be seen that the discharge inpeadance during the bulk of the energy deposition varies 
between 5 D and 10 D. This variation suggests an optima value of -7 D. for the circuit 
inpeadance. It should be noted that the discharge iipeadance will vary with gas conposition and 
optinisation of the discharge circuit should be undertaken after the final gas nixture has been 
decided. 

4.2.3 Summary 

Voltage and current pulse neasurenents of the LP-140 discharge were obtained and fron these were 
derived the electrical energy deposited into the gas, the electrical to optical conversion 
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Figure 4.6. The Discharge Inpeadance. 

efficiency, the input energy dependence on E/N and the discharge inpeadance. It was found that 
3.4 J was deposited into each discharge, which had an inpeadance of (5-10) Cl A poor electrical to 
optical conversion efficiency of -4 \ was found to be due to ringing of the voltage pulse which 
resulted in the electrical energy being deposited at too high a value of E/N for efficient excita- 
tion of the relevant COa and Na vibrational levels. 

4.3 Power supply considerations 

4.3.1 Introduction 

This section provides an outline of assistance provided to Philip Thistlethwaite, a GSRP student who 
worked with the LP-140 during the later half of 1990. His task was to design a new high-voltage 
power supply for the LP-140. 

4.3.2 Design Criteria 

The following Hathcad docunent was provided to Philip as an initial specification requirenent. A 
noninal design pulse rate of 20 Hz was chosen to enable the laser to be run at high pulse rates if 
required. The pulse energy was specified at 200 J and voltages up to 40 kV to enable higher gas 
pressures to be used in the laser head. This would result in shorter pulses than the (5-15) (.ts i 
pulses currently obtained, together with a larger output energy. As the design of high voltage power 
supplies becones nore difficult with increasing voltage, the naxinui voltage required was set at 
20 kV but with the request that 40 kV would be desirable if he found it possible. 
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Power Supply Design Criteria 

This determines some of the basic requirements of the high voltage power 
supply based on the required inputs to the laser which are;- 


Maximum 

voltage 

(V) 

V 

:= 20000 

Maximum 

pulse 

energy 

(J) 

max 

E 

:= 200 

Maximum 

pulse 

rate 

(HZ) 

max 
PRF : 

= 20 

Average 

power 


(W) 

P 

: = PRF • E 





av 

P 

av 

max 

3 

= 4 10 


Allow a 100/nS period for discharge and thyratron recovery, during which the 
charging unit is off, therefore the time available to pulse charge a capacitor 
is given by 

-4 

t := 1- 10 S 

discharge 

1 

t := t 

charge PRF discharge 


t = 0.0499 S 

charge 

Thus allow a 49mS charge time 


t := 0.049 S 

charge 

This gives an average charging current for the chosen pulse energy of 

E 

max 

I := 

av V • t 

max charge 


giving I = 0.204 A 

av 
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V * 


The capacitor will charge exponentially where the charging current will be 
given by 


1=1 exp(-ct) 
pk 


where I is the peak charging current, c is the time constant and t is the 
pk 

time. 

When the capacitor is fully charged the current will fall to zero, however due 
to the exponential decay the point at which this occurs can be said to vary, 
depending on how accurately the capacitor is to be charged. If we assume that 
we require the capacitor to be fully charged to better than 0.5% then this 
will give a time and time constant product of:- 

ct := -ln(0.005) 

ct = 5.298 

This gives for the time constant, c assuming the charging time t above. 

charge 


ct 


t 

charge 


-1 

C = 108.129 S 


We now need to calculate the peak charging current. The average current 
obtained above can be regarded as the integration of the charging current over 
the charging time, divided by the charging time ie 


I 

pk 


I • t 
av charge 


t 

charge 

exp(-c-t) dt 


J 0 


I = 1.087 A 
pk 
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,7e can then plot the charging cycle 


t := 0,0.001 . .0.05 

I(t) := I -exp(-c-t) 
pk 



Time elapsed since start of 
charging cycle (S) 
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4.3.3 Circuit Outline 

Figure 4.7 shows a suggested circuit outline for the power supply which was provided to Philip. The 
input terminals on the left-hand side assume a dc source, Vi„ which charges up capacitor Cl. 
Capacitor C2 is resonantly charged through inductor LI to a voltage 2xVi„. To ensure negligable 
ripple on C2, C1>C2. When C2 is fully charged thyristor, T1 is triggered to discharge C2 through the 
high-voltage step-up transformer, TR1. Due to the large size of Tl, it requires a finite time to 
turn on and saturating inductor SL1 limits the current through the thyristor until turn on has been 
completed, at which point SL1 saturates to provide negligible resistance to the discharge pulse. 



Figure 4.7. An Outline PSD Circuit. 

The transformer, TR1 is shown as a distributed core transformer. This enables a large step-up ratio 
to be obtained whilst maintaining a low inductance for the transformer. It is possible that a single 
core high voltage transformer could be used. After being stepped up the pulse appears on capacitor 
C3. Diode, D1 prevents reverse current fluctuations from damaging the thyristor, which is vulnerable 
to such pulses. By choosing a suitable value of C3 and charging and discharging C2 rapidly an 
apparently d.c. high voltage is obtained at the output terminals. The resistive divider on the 
output enables the output voltage to be monitored. This monitor feeds back to the thyristor switch- 
ing circuit (not shown). For a given load, varying the thyristor pulse rate will result in a 
variation of the output voltage. 
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4.3.4 Pulse Trans f oner Design 

The Kathcad docuient following the references enables a design for a single core low inductance 
pulse transfoner to be evaluated. This docuient was prepared following a request for a transformer 
design froi the student. The values in the docuient are those for a transfoner designed previously 
by the author' 4-2 ’. These values were used to test the correctness of the calculations. As no 
specifications for the transfoner were recieved from the student a final design was not completed. 
The theory for the design of low inductance high-voltage transfoners is soievbat complex and long 
and the reader is therefore referred to the standard texts by Snelling' 4-3 ' 4-4 ’ and the paper by 
Baker' 4-5 ’ which provide a coiplete review of the subject. The sunary of data at the end of the 
docuient enables the optiiui transfoner to be chosen. 

4.4 References 

4.1) Lowke, J.J., A.V. Phelps i B.W. Irwin; "Predicted electron transport coefficients and operating 
characteristics of COa-Ha-He mixtures"; J. Appl. Pbys.; 44; 4664 (1973) 

4.2) Sylvan, A., G.D. Spiers, A. Johnstone & LG. Harrison;"Coipact high repetition rate nagnetically 
switched transverse, electric, atmospheric CCh laser"; Rev. Sci. Instrum.; 61; 1821 (1990) 

4.3) Snelling, E.C.;"Soft Ferrites: Properties and Applications"; Butterworths (1988) 

4.4) Snelling, E.C. i A.D. Giles; "Ferrites for Inductors and Transfoners"; Research Studies Press Ltd. 
(1983) 

4.5) Baker, H.J., P.A. Ellsmore i E.C. Sille, J. Phys. E: Sci. Instrum.; 21; 218 (1988) 
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1000 


V 


Input Parameters 
Peak Input Voltage 
Peak Output Voltage 

Voltage Ratio 


Pulse Transformer Design 
Gary D. Spiers 


VI 

pk 

VO := 30000 V 

pk 

VO 

pk 

Vratio := 

VI 

pk 


Pulse Width 
Pulse Energy 

Assume pairs of 3C8 FX3845 ferrite E-cores 
Assuming the cores are biased then the maximum 
available flux swing is 

Effective area of magnetic path 

Effective magnetic path length 

Effective magnetic volume 

Available winding height 


-6 

r := 8-10 
E := 4.32 


(SB 

A 

e 

1 

e 

v 

e 

Wh 


:= 0 . 85 T 

-6 2 

: = 532-10 m 

:= 0.147 m 

-5 3 

- = 7.82-10 m 
-3 

:= 11-10 m 


If we stack m cores together and vary nl , the number of turns on the primary 
between l and 7 then the condition to avoid saturation is given by:- 1 


nl := 1 . . 7 


VI ■ T 
pk 

m := 

nl <SB-A ■ nl 

e 

Make allowances for need to have an integral number of cores 

m : = m + l - mod m , l' 
nl nl [ nl 

and the number of secondary turns is given by 

n2 := nl- Vratio 
nl 
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M 



We require an integral number of turns on the secondary 


n2 := n2 - mod 
nl nl 


n2 , 1 
nl 


Now calculate the primary inductance 
Relative Pulse Permeability of Core 

Permeability of Free Space 


:= 720 


-7 


:= 4 • 7T • 10 


H/m 


Lp 


nl 


/j • /j -A • m • nl 
0 p e nl 6 

10 /jH 



-4 


t 

:= 3.8-10 

m 

1 

-5 


t 

o 
< — 1 

in 

♦ 

II 

m 

2 

-5 


c 

:= 2.4-10 

m 


-2 


w 

:= 1.9- 10 

m 


Next the leakage inductance assuming a tape wound core 
Thickness of material in primary winding 

Thickness of material in secondary winding 

Thickness of interwinding insulation 
Winding Width 

The thickness of the primary winding is thus given by: — 

tp := ft nl + c- (nl + 1)] m 
nl L 1 

and the mean length per turn on the primary is given by:— 

-2 -2 

lp := 2.74-10 ■ 2 • m - + 2-10 -2 + 7 T-tp m 

nl nl nl 

Similarly for the secondary the thickness of the winding is given by:- 


ts 


nl 


~t • n2 

+ c- 

"n2 + l' 


2 

nl 

nl 



m 


and the mean length per turn is: 


-2 -2 

Is := 2.74-10 ■ 2 ■ m + 2-10 -2 + 4 • tp + n-ts m 

nl nl nl nl 
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thus providing an overall mean length per turn of:- 


nl • lp + n2 ■ Is 

nl nl nl 


nl 


nl + n2 


nl 


and a total winding thickness, tw 


nl 


tp + ts 
nl nl 


10 mm 


The leakage inductance is given by:- 


-1 2 nl 

11 := 4- 7T- 10 • nl 

nl w 


t • nl + t • n2 
1 2 nl 


- + c- fnl + n2 

L nl 


MH 


From the pulse energy and the input voltage we can calculate the input 
capacitor size H 


Input Capacitor 


2- E- 10 

Cl := nF 


Cl = 8.64 10 


nF 


VI 


pk 


If we assume the inductance of the circuit in series with the primary and 
secondary coils is of the same magnitude as the leakage inductance and also 
assume that the capacitor on the output side of the transformer is iust 
sufficient to hold the energy transfered across the transformer then we can 
calculate the transfer times for the various configurations. " 


9 

2- E- 10 

Output Capacitor C2 := nF 

Referred to primary 2 

VI 

pk 


Series inductance 


3 

C2 = 8.64 10 nF 


Ls := 11 /iH 

nl nl 


Transfer Time 


t 

nl 


TT • I ["11 

'i L nl 


+ Ls 

nl 



C1-C2 

. j 

Cl + C2 


-3 


• 10 


MS 
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The actual frequency is thus given by 


n 

oj := rad/s and 

nl -6 

t • 10 
nl 

and the actual flux swing in each 


1 

f := 

nl -6 

2 * t -10 
nl 

core is given by 


Hz 


VI 

pk 

dB := T 

nl nl-ra -A -w 

nl e nl 


The energy losses in the transformer can be divided into two distinct 
one due to the winding and the other due to the magnetic material. 

Magnetic losses 

There are two sources, hysteresis loss and eddy current loss. 
Hysteresis Loss 


Steinmetz Exponent 

Proportionality Constant 

Energy loss per pulse due to hysteresis 


n := 2.5 

-2 

k := 2.154- 10 


Eh 

nl 


m • V 
nl e 


• k- 



n 6 
■ 10 


mJ 


f> := 2 ftm 

f 

R := 16 


2 2 
n • m -A • dB -V 

nl e nl e 3 

Ee : 1 0 mJ 

nl f> ■ 6- r 

f 


Energy loss per pulse due to eddy currents 
Resistivity of core material 
Filling factor for core geometry 


sources 


26 



Winding Losses 

Resistivity of copper winding 


-8 


1° := 1.694- 10 flm 

c 

The resistances of the primary and secondary coils are thus:- 


Rldc 

nl 


p -nl-lp 
c nl 

n 

t • w 
1 


and 


R2dc 

nl 


p • n2 -Is 
c nl nl 
n 

t • w 
2 


respectively. The peak currents in the primary and secondary can be 
estimated by:- 


IP 

nl 


and 


VI 

pk 


to • fll 
' nl [ nl 


-6 


+ Ls 


nl 


10 


A 


IS 

nl 


VO 


pk 


co • rn 

nl |_ nl 


+ Ls 

nl 


-6 2 
• 10 • Vratio 


A 


respectively. The energy loss per pulse in the core windings due to the d.c 
resistance can then be calculated 


Eres 

nl 


2 

IP • Ride 
nl nl 


2 

+ IS -R2dc 
nl nl 



mJ 


In general the a.c. resistance is higher than the d.c. resistance at high 
frequencies due to eddy currents in the windings. These manifest themselves as 
the skin effect and the proximity effect. 
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Skin effect 



For these merit fuctions less than 2 the skin effect can be ignored. 
Proximity Effect 

For tape wound cores the energy loss per pulse is given by:- 

2 2 
TT • dB • w 

nl I" 3 3] 3 

Epe := nl-lp -t + n2 -Is t -10 mJ 

nl 24 • t ■ f> L nl 1 nl nl 2 

c 

The total energy loss is given by 

Etot := Eh + Ee + Eres + Epe mJ 

nl nl nl nl nl 


The results from these calculations are tabulated and summarised in 
graphical form on the following pages. 
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A summary of the variables in the following tables and their respective units 
are : - 


nl , the number of primary turns (no units). This can be considered the primary 
variable, all other variables having been derived with respect to nl. In the 
following variable summary the nl subscript has been dropped to reduce the 
amount of space required. 

m, the number of core pairs required (no units). 
n2, the number of secondary turns (no units). 

Lp, the primary inductance (uH). 

11, the leakage inductance (uH). 

t, the pulse transfer time (us). 
dB, the calculated flux swing in the core (T). 

IP, the peak current in the primary (A). 

IS, the peak current in the secondary (A), 
tw, the total winding thickness (mm). 

Ml, the skin effect merit fuction for the primary (no units). 

M2, the skin effect merit function for the secondary (no units). 

Eh, the energy loss due to hysteresis loss in the core (mJ). 

Ee, the energy loss due to eddy current loss in the core (mJ) . 

Eres, the energy loss due to the d.c. resistance of the windings (mJ) . 

Epe, the energy loss in the windings due to the proximity effect (mJ) . 

Etot, the total energy loss in the transformer (mJ) . 



m 

n2 

Lp 

n 

nl 

nl 

nl 

nl 

nl 

1 

18 


30 


58.9 


0.08 

2 

9 


60 


117.9 


0.32 

3 

6 


90 


176.8 


0.77 

! 4 

5 


120 


262 


1.61 

5 

4 


150 


327.4 


2.68 

: 6 l 

3 


180 


353.6 


3.86 

7 i! 

3 ; 

210 


481.3 

! | 

6.28 


dB IP IS 

nl nl nl nl 


2.6 


0.08 


5312.3 


177.1 

5.3 


0.17 


2584.1 


86.1 

8.1 


0.27 


1669.6 


55.7 

11.7 


0.35 


1159.9 


38.7 

15.1 


. 0.45 


897 


29.9 

18.1 


0.6 


748.2 


24 . 9 

23 . 1 


0.66 


586.4 


19.5 


tw 



Ml 


nl 

nl 

1 


2.567 

2 


1.79 

3 


1.439 

; 4 


1.199 

; 5 

! 

1.055 

i 6 


0.963 

! 7 


0.853 


M2 

nl 

Eh 

nl 

0.169 


0.1 

0.118 


0.2 

0.095 


0.4 

0.079 


0.6 

0.069 


0.9 

0.063 


1.4 

0.056 


1.8 


Ee Eres 
nl nl 



262 

.4 

133 

.7 

90. 

8 

72. 

5 

59 

47. 

8 

44. 

3 


Epe 

nl 


23.6 

104.2 

260.5 

499.4 

868.8 

1477.9 

2075 


1 


Etot 

nl 

287 . 3 1 

239.4 ! 

~ 353 i 

574.1 

930.4 
1528.8 

2123.1 
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All graphs are a function of nl , the number of primary turns. 


20 

Number of 
Core Pairs m 

nl 


1 

1 nl 7 



Number of 
Secondary 
Turns 



0.7 

Flux 

Swing dB 
(T) nl 

0 

1 nl 7 



Primary 

Inductance 

(uH) 



Leakage 

Inductance 

(uH) 



Transfer 

Time 

(US) 


30 


t 

nl 
0 

1 nl 7 



Primary Peak 
Current 
A 


6000 

IP 

nl 

0 


K 

j ! 


1 r * 


nl 


Secondary Peak 
Current 
A 



Winding 

Height 

mm 



30 









Core 

Hysteresis 
Loss ( mJ ) 


Eh 



Resistance 

Loss 

(mJ) 



3000 

Total Energy 

is: 

iJ 


Core Eddy 
Current Loss 
(mJ) 


Proximity 
Effect Loss 
(mJ) 






5 Frequency Chirp 

5 . 1 Introduction 

There are two contributions to frequency chirp during the laser pulse, one due to the decaying electron 
density' 5 ' 11 in the discharge and the other due to theraal changes in the laser gas' 5 ' 21 . The theraal 
effect or Laser Induced Hediun Perturbation (LIMP) is the sain concern. 

5.2 Analysis 

Frequency data is collected within the laboratory on a HP tiie and frequency analyser and stored on an 
IBM coipatible coiputer. The data is analysed using the first of the two Hathcad documents which follow. 
The docuient is self explanatory. The saiple file analysed in the document is a file obtained froi the 
LP-140. It can be seen that the total pulse length is 13 jjs long with a total frequency chirp of 
2.7 MHz consisting of an initial downsweep due to the decaying electron density and then an upsweep due 
to the LIMP contribution. Within the center of the pulse there is a 5 ^is portion with a frequency chirp 
of 230 kHz and this has been plotted in the second window. This docuient has been nade available to Dr. 
Jaenisch to enable him to analyse the chirp results froi the LP-140. 

5.3 Modeling 

The LIMP within the laser has been nodeled in Hathcad using the theory of Willetts et al.' 5 ' 21 . Froi the 
second Hathcad docuient we see that for a pulse length of 13 ps and a pulse energy of 1 J a frequency 
chirp of 2.5 MHz is obtained. As the LIMP theory is only an approximation, the degree of agreeient 
between the neasured chirp and the calculated chirp is very good, with a less than 101 difference 
between the two values. 

5.4 References 

5.1) Willetts, D.V. & H.R. Harris;"A Plasia Effect in Injection-Mode Selection of TEA C0 2 Lasers"; Appl. 
Phys. B; 33; 91 (1984) 

5.2) Willetts, D.V. & H.R. Harris; "Scaling Laws for the Intrapulse Frequency Stability of an Injection 
Mode Selected TEA CO 2 Laser"; IEEE J. Quant. Electron.; QE-19; 810 (1983) 
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H.P. DIGITISER CHIRP DATA PLOTTING AND ANALYSIS 


Gary D. Spiers 

The heterodyne beat frequency obtained to measure the frequency chirp in the 
LP-140 is digitised by a HP digitiser and the data passed to a PC where it is 
stored for later use. Each file consists of 201 data points where each point 
consists of the point number, chirp frequency and time interval since the 
previous measurement respectively. 

This document extracts the time and frequency information from the data file 
and plots two windows, the first is the total chirp over the entire pulse 
whilst the second is a user defined sub-window of the first. The frequency 
change and time duration of each window is also calculated. 


i := 0 ..200-3 - 1 This is the index for reading in the data. 


data 

i 


READ 


chirp 

dat 


This reads the data. 


j := 0 ..199 Set up index for extracting required information. 


The following two statements extract the chirp and time data respectively. As 
the original data has units of Hz and s respectively the data is also 
converted into kHz and ms respectively. 


data 

j- 3 + 1 

chirp := 

j 3 

1 - 10 


time 

j 


data 

j- 3 + 2 


-6 

1 - 10 


Now the two windows are set up. As the digitisation period is much longer than 
the pulse length, the first window is used to display the portion of relevance 
and the second to highlight particular areas of interest. The variables 
cstart, cend, sstart and send define the size of each of the windows, the 
units being the number of points . 


cstart := 22 cend := 77 sstart := 31 send := 51 


The following statements set up the index range for each window and copy the 
relevant data points from the original data arrays. 


k := 0 ..(cend - cstart) 

c := chirp 
k k+cstart 

d := time 
k k+cstart 


1 := 0 ..(send - sstart) 

e := chirp 
1 1+sstart 

f := time 
1 1+sstart 


The following statements find the maximum and minimum frequency in each window 
and calculate the total frequency chirp and time duration of each window. The 
frequency is also offset to provide a minimum frequency of 0 kHz to enable 
convenient plotting. 

mine := min(c) mine := min(e) 

maxc := max(c) maxe := max(e) 


(Schirpl := maxc - mine 


6chirp2 := maxe - mine 


c := c - mine 
k k 


timel 



d 

k 



Finally plot the data and display the calculated data. 


3000 


k 


0 


0 

k 

55 


timel =13 

AlS 


(Schirpl = 2700 

kHz 
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Laser Induced Medium Perturbation 
Gary Spiers 


Gladstone-Dale coefficient, 

Equivalent frequency of lower laser level, 
41.9 THz for 10 jum band, 38.7 THz for 9 /im 


-4 

K := 2.11-10 

12 

fO := 41.9-10 


Output pulse energy, (J) E 
Molar gas constant, r 
Mode spot radius, (M) a 
Cavity length, (M) L 
Specific heat capacity at constant volume, Cv 
Optical Pulse length, (/iS) r 
Time,t (liS) t 


1.0 


8.314 

-2 

1.25-10 


2.5 


16.97 

13 

0,0.1 . . r 


f(t) 


3 

2 • K- f0-E-R-t 


4 15 

3 • n- a • L- Cv- r • 10 


The factor of 10 A 15 
accounts for the 
conversion of t and t from 
MS to S and Hz to kHz 


3000 


Frequency 
Chirp f(t) 
(kHz) 


0 

o t 13 

Time,t /iS 



Total frequency chirp, f(r) = 2545 kHz 
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6 Laser Linewidths 

6.1 Introduction 

The laser linewidth is inportant in determining the likelihood of obtaining single axial node operation. 
For pressures above 20 torr the COz linewidth is a conbination of two lineshapes. The first is due to 
Doppler broadening, which is the dominant nechanisi at low pressure, and has a Lorentzian lineshape. The 
second and usually uuch larger contribution arises fron collision broadening which has a Gaussian 
profile. The combined profile is known as a Voigt profile. For pressures above -50 torr however the 
collision broadening is so uuch larger than the Doppler broadening that the Doppler can be ignored for 
most situations. The Hathcad docuient that follows calculates a Gaussian lineshape to approxiiate the 
Voigt profile. This enables a conparison between the LP-140 and the lasers proposed by the LAWS contrac- 
tors to be nade. For the LAWS lasers this approximation is inconsequential, however for the LP-140, 
which operates at a nuch lower pressure, there will be a slight error. 

6.2 LP-140, LMSC/AVCO and GE/STI Laser Comparison 

Table 6.1 lists a conparison of the results obtained for each of the lasers. LP140a is the LP-140 with a 
short pulse gas mixture and LP-140b with a long pulse gas nixture. 


Laser 

Cavity 

P 

C0 2 

N 2 

Linewidth 

Axial Mode 

Ho. modes 


Length (n) 

(torr) 

(?) 

(?) 

FWHH (HEz) 

Sepn. (MHz) 

under FWHH 

LP-140a 

2.38 

45 

52 

29 

343 

63 

5.4 

LP-140b 

2.38 

45 

13 

15 

291 

63 

4.6 

LMSC/AVCO 

2.20 

380 

25 

25 

2215 

68 

32.5 

GE/STI 

3.00 

380 

17 

33 

2151 

50 
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Table 6.1. LP-140, LMSC/AVCO and GE/STI Laser Mode Conparison 


From the table it can be seen that there is a distinct difference between the LP-140 and the LAWS 
lasers. This is due to the difference in operating pressures, the LAWS laser linewidths are dominated by 
the collision broadening. The number of axial modes under the FWHH of the linewidth indicates the 
likelihood of axial mode hopping, the greater the number the more likely mode hopping will occur. Thus 
in principle the LP-140 should easily hold a single axial mode, however due to its construction it 
exhibits considerable thermal drift which negates the narrow linewidth advantage. A cavity length 
increase of 5 nm is sufficient to change the axial node yet perspex and aluminium, the two main compo- 
nents of the LP-140 construction have thermal expansion coefficients of 85 (.imK -1 and 23 unK' 1 respec- 
tively. For a mechanical resonator length of -l n this implies a temperature stability of -(0.25-0.06)K 
is required to keep the resonator stable. 

6.3 The Mathcad Document 

As mentioned previously this document only approximates a Voigt profile, however modifications are in 
hand to calculate the exact profile. 
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Laser 


This calculation normalises 


Linewidth Calculation (LP-140) 

Gary D. Spiers 

the gain shapes to one. 


Cavity length, (m) L 

Speed of light in vacuo, (m/s) c 

Wavelength of interest, (m) X 

c 

Rotational line number J 

Total gas pressure, (torr) P 

Fractional C02 concentration C02 

Fractional N2 concentration N2 

Fractional He concentration He 

Gas Temperature, (K) T 

Avogadro's number, (/mol) NA 

Boltzmann's constant , (J/K) k 

Molecular mass of C02, (Kg/mol) MC02 

No. of points calculated imax 


Doppler Broadening (FWHM) 


<5vd 


8- In ( 2 ) -k-T 


c 


j r MC02] 2 

i — i - c 

■I na J 


x 

c 


7 

<Svd = 5.289 10 Hz 


Collision Broadening (FWHM) 


6 vc 


<S vc 


20.660- (7.520 - 


sjT 

8 

= 2.903 10 HZ 


0 . 059 ■ J) • P 

( C02 + 0.73-N2 +0.64 
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:= 2.38 

:= 2.99792458- 10 
-6 

:= 10 . 6-10 

:= 20 

:= 45 

:= .5193 

:= 0.2909 

:= 1 - C02 - N2 

:= 300 

2 

:= 6.0221367- 10 
-2 

:= 1.380658-10 
-3 

:= 44.0098-10 
:= 400 


6 

He) • 10 


Total broadening ( FWHM) 
6vc := <Svc + 6vd 
8 

6 vc = 3.432 10 Hz 


Axial Mode Seperation 


<5vax : = 


2 • L 


<5vax = 


nax : = 


6.298 10 
6 vc 


Svax 


Hz 


Lineshape calculation 


i : = 1,2 . . imax 


:= 282813 + — 

L 20 


c 

v 


6vc 


ac : = 
i 



c 


V 



i 

X 



L c J 



+ (Svc 


8 

• 10 
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Axial mode plotting calculations 


1 


imax + l 


int := i - mod (i, step) 
i 


d := if Tint ^ int ,0.5,0 
i i i+1 


j = = 

int := 0 

j 

imax- 6vax 

step := — 

v - v 

imax 1 


-6 

Svc := <Svc- 10 
(Conversion to MHz) 


LP-140 Linewidth 



Frequency 



Minimum 


Maximum 



13 

13 

Frequency, (Hz) 

v = 2.82813 

10 

V = 2.82833 10 


1 


imax 



-5 

-5 

Wave length, ( m ) 

X = 1.05996 

10 

X = 1.06004 10 


imax 


1 


Linewidth (FWHM), <5vc = 343.2 MHz 

No. of axial modes under the FWHM, nax = 5.4 
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7 Appendix ( 1 ) 


PROGRAM ANLYSVI 
Gary D. Spiers 

Center for Applied Optics, University of Alabama, 

Huntsville, Al. 35899, (205) 895 6030 

This prograi processes digitised voltage and current pulses obtained 
froi the LeCroy 9450 350 MHz digitising oscilloscope. The data has been 
previously saved to disc using prograi LECROY, written by Phil Kroiis of 
Boeing Coiputer Services. This prograi stores the data in two fons, the 
raw data direct froi the digitiser (contained in the directory c:\lecroy\ 
lcrdata) and decoded data (in the directory c:\lecroy\lcpdata). This 
prograi uses the processed data as input. The data in these processed files 
is stored as follows 


Vertical axis gain per digitisation point (Volts or Amps) 

Vertical offset (Volts or Amps) 

Timebase interval between digitisation points (s) 

Tiiebase offset with respect to the oscilloscope trigger pulse (s) 
Timebase value (s/div) 

(A blank line which may eventually contain info, added at storage time) 

!! ri f? n n it ft ti n n n 


ri 


I! 


!1 


n 


!1 


It 


1st vertical value 


final vertical value 
end of file 


The output files (except for the energy dependence on E/H) are stored in a 
siiilar fomat:- 

Time between digitisation points (s) 

Time offset froi trigger point to first digitisation point (s) 

1st calculated vertical value in appropriate units 


final calculated value. 
End of file. 


' The energy as a function of E/N file consists of 11 values of energy, one 
' for each E/H subdivision. 

' START OF PROGRAM 

' Dimension the filename arrays and the array for the energy as a function 
' of the reduced electric field data. 
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DIM filel AS STMHG * 40 
DM file2 AS STMHG * 40 
DM file3 AS STMHG * 40 
DM file4 AS STMHG * 40 
DM file5 AS STRMG * 40 
DM file6 AS STMHG * 40 
DM file7 AS STMHG * 40 
DM files AS STMHG * 40 
DIH file9 AS STMHG * 40 
DIH fileio AS STMHG * 40 
DIH filnai AS STMHG * 15 
DIH ebin(l TO 11) AS DOUBLE 

' Declare the variables used in the program 

DIH dl, d2, d3, d4, d5 AS SIHGLE 
DIH vqain, voffset, hstep, hoffset, timebasl AS DOUBLE 
DIH igain, ioffset, tstep, toff set, tinebas2 AS DOUBLE 
DIH vpoint, ipoint, efn, nrgstep, oldnrg, energy AS DOUBLE 
DIH iiipead, efnval, tempry, edrift, electron AS DOUBLE 
DM count AS IHTEGER 


' Assign values to constants 


C0HST pressure = 50 
COHST temp = 100 
COHST R = 8.315 
COHST elecsep = 4 
COHST elecwid = 4 
COHST e = 1.602177D-19 

'Assign filenames for input and 


'Total gas pressure in torr 
'Gas temperature in centigrade 
'Universal gas constant 
'Discharge gap in ci 
'Discharge width in ci 
'Charge on an electron 

of data 


filel = "c:\lecroy\lcpdata\" 
file2 = "c:\lecroy\lcpdata\" 
file3 = "c:\lecroy\lcpcalc\energy.dat" 
file4 = "c:\lecroy\lcpcalc\iBpead.dat" 
file5 = "c:\lecroy\lcpcalc\efn.dat" 
file6 = "c:\lecroy\lcpcalc\edrift.dat" 
file7 = "c:\lecroy\lcpcalc\electron.dat" 
files = "c:\lecroy\lcpcalc\voltage.dat" 
file9 = "c:\lecroy\lcpcalc\current.dat" 
filelO = "c:\lecroy\lcpcalc\nrgven.dat" 

'Clear the screen and obtain the names of the digitised voltage and 
'current pulses to be analysed. 


SCREEH 9 
CLS 

IHPUT "Voltage filename "? filnam 
filel = LEFT$( filel, 18) + filnam 
IHPUT "Current filename "? filnam 
file2 = LEFT$(file2, 18) + filnam 


'Open up the files for input and output respectively 
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'Input voltage 
'Input current 


OPEN filel FOB IHPOT AS #1 
OPEN file2 FOE IHPOT AS J2 
OPEN file3 FOE OOTPOT AS #3 
OPEN file4 FOE OOTPOT AS ft 
OPEN file5 FOE OOTPOT AS /5 
OPEN file6 FOE OOTPOT AS J6 
OPEN file7 FOE OOTPOT AS t 7 
OPEN file8 FOE OOTPOT AS it 
OPEN file9 FOE OOTPOT AS /9 

' Bead in the voltage gain and offset, digitisation tine interval, pre/post 
'trigger delay interval and oscilloscope tiiebase setting for the voltage 
'pulse followed by five blanks to separate the header fro* the data 

INPOT tl, vgain, voffset, bstep, hoffset, ti*ebasl 
INPOT tl, dl, d2, d3, d4, d5 

'Eepeat the above procedure for the current pulse to obtain current gain and 
'offset, digitisation ti*e interval, pre/post trigger delay interval and 
'oscilloscope tinebase setting for the current pulse. 

INPOT 12, igain, ioffset, tstep, toffset, ti*ebas2 
INPOT tl, dl, d2, d3, d4, d5 

' As the two pulses were digitised simultaneously the digitisation intervals 
' and offsets should be the sa*e. If there is a difference this implies an 
' error in the data or possible incorrect pairing of the voltage and 
' current data files. If they are not the saie then issue a warning and stop 
'the progra*. 

IF toffset <> hoffset OE hstep <> tstep THEN 
PBINT "Voltage and Current files are incompatible !" 

CLOSE 
STOP 
END IF 

' As all of the output files except the energy deposition as a function of the 
' reduced electric field are tine dependent, write the digitisation interval 
' and trigger offset to the top of each file (except Energy as fn. of E/N). 

WEITE #3, tstep 
WEITE #3, toffset 
WEITE #4, tstep 
WEITE #4, toffset 
WEITE #5, tstep 
WEITE t 5, toffset 
WEITE t 6, tstep 
WEITE t 6, toffset 
WEITE tl, tstep 
WEITE tl, toffset 
WEITE if 8, tstep 
WEITE tt, toffset 
WEITE f9, tstep 
WEITE #9, toffset 
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' Set up a loop with counter to increientally read in tbe data points and 
' calculate the required intonation. 

count = 0 
CIS 

LOCATE 10, 28 

PRIST "Processing point no." 

DO 

LOCATE 10, 49 
PRIST count 
INPUT #1, vpoint 
IRPUT #2, ipoint 

'Firstly take the raw data points and convert then to units of volts and anps 

vpoint = -(vpoint * vgain - voffset) 
ipoint = (ipoint * iqain - ioffset) 

' Then store these voltage and current values in their output data files so 
' that they nay be plotted 

WRITE #8, vpoint 
WRITE #9, ipoint 

' Now calculate the reduced electric field, E/N and store the value 

efn = ((273.15 + tenp) * R * vpoint) / (80.2853 * 10 A 18 * pressure * elecsep) 
WRITE #5, efn 

' Multiplying the voltage by the current and tinestep provides the input 
' energy for each discrete digitisation interval. This energy per step is 
' stored in an array for use later to evaluate the energy as a fn. of E/B. 

' The individual packets of energy are sunned to provide the total energy 
' input into the discharge and stored to file, this provides a record of 
' the tine evolution of the energy into the discharge. 

nrgstp = vpoint * ipoint * tstep 
IF count = 0 THEM 
energy = nrgstp 
ELSE 

energy = oldnrg + nrgstp 
END IF 

oldnrg = energy 
WRITE /3, energy 

' The discharge inpeadance at each digitisation point is found by dividing 
' the voltage by the current. However in the initial and final parts of the 
' digitisation both the current pulse is zero and the division results in 
' an infinite inpeadance which creates an overflow error. Additionally the 
' current fluctuates by +/- 1 digitisation point (an artifact of the 
' digitisation process) about the zero level to provide spurious values 
' of the inpeadance which are artificial and due entirely to the digitisation. 

' To elininate this traps for the nagnitude of the discharge current <0.1 A 
' (c.f. peak values of -500+ A) have been set and when this occurs the 
' iipeadance is set to 200 ohns which is sufficiently large with respect to 


43 



' the actual impeadance after breakdown has occured to be representative of 
' the actual infinite value prior to breakdown. The resulting iipeadance 
' value is saved to file. 

IF ABS(ipoint) < .1 THEN 
iipead = 200 
ELSE 

iipead = vpoint / ipoint 
END IF 

IF iipead <= 0 THEN iipead = 200 
WRITE J4, iipead 

' This routine subdivides the input energy as a function of the E/N value. 

' The E/N is divided into 10 boxes, each representing a step of lxl0 A (-16) 

' Vci A (3) in the range (0->10)xl0 A (-16) Vci A (3) plus an eleventh box for 
' all values greater than 10xlO A (-16) Vci(3). As only the final total for 
' each box is required no saving to file is done at this point. 

efnval = ABS(efn) * 10 A 16 
IF efnval < 10 THEN 

ebin(INT( efnval + 1)) = ebin(INT( efnval + 1)) + nrgstp 
ELSE 

ebin(ll) = ebin(ll) + nrgstp 
END IF 

' The following calculates the electron drift velocity and electron density 
' in the discharge. The electron density depends on the electron drift which 
' depends in a non-linear fashion on the reduced electric field E/N. The 
' paper "Predicted electron transport coefficients and operating 
' characteristics of C02-N2-He laser mixtures" by J.J. Lowke, A.V. Phelps and 
' B.W. Irwin in J. Appl. Phys., 44, 10, 4664-4671, (1973) calculated this 
' dependence for several gas mixtures. From the results presented we can see 
' that :- 

r 

' log(drift velocity) = b + k.log(E/N) 

/ 

' approximately. Allthough crude this approximation enables an order of 
' magnitude estimate of the electron drift velocity to be made simply. This 

' technique has been used within this program with b=17.044 and k=0.658. 

/ 

' A trap is used to prevent a value of zero volts resulting in an overflow 
' error. Finally the calculated data is saved to file. 

IF vpoint <> 0 THEN 
tempry = LOG(ABS(efn)) / L0G(10) 
edrift = 10 A (17.044 + .658 * tempry) 
electron = (ipoint / (elecsep * elecwid * e * edrift)) 

ELSE 

electron = 1 
edrift = 1 
END IF 

WRITE /6, edrift 
WRITE il, electron 
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' Increment the loop counter and if either of the input data files is 
' eipty then leave the loop and close all the files. 

count = count + 1 
LOOP UNTIL EOF(l) OR E0F(2) 

CLOSE 

' Finally the output energy dependence on E/N data mist be written to file 
' and then end the program 

OPEN filelO FOR OUTPUT AS #10 
count = 1 

DO 

WRITE #10, ebin(count) 
count = count + 1 
LOOP UNTIL count = 12 

CLOSE 

CLS 

LOCATE 10, 28 

PRINT "Analysis successfully completed. n 
END 


45 


8 Appendix (2) 

Idealised Laser Mode Profiles 


This program oalovil^es the Rayleigh range for. a given Gaussian beam 

together with the spot size and mode profile at any given location. The 

calculated is the idealised rectangular mode profile of the beam at 

1 ° C ; tl ? n for values of 1 and ». TEM indices, up to 3 
select the initial parameters:- H ° firstly 


Minimum spot size 
-3 

w := 3 • 10 
0 


Wavelength 

-6 

X := 10.6-10 


Refractive Index 
n := l 


Calculate the Rayleigh range at the beam waist, given by:- 


2 n 

z : = tt • w - 

R ox 


Z = 2.667 
R 


Define the beam waist location as position z=o then 
calculation of the spot size and mode profile :- 


chose 


a position for the 


z := 0 
0 


z := 5 


Calculate the spot size at the chosen location using: - 



2 

2 

2 

w 

i + 

0 

2 


z 


R 


w = 0.006 


Then calculate the radius 
using :- 


of curvature of the beam at the chosen location 


R 


z ■ 


1 + 



R = 5 
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Now start on the mode profile by chosing values 
indices, and the electric field, E. The mode is 
of points. 


for 1 and m, the mode 
calculated using a 50x50 arra' 


1 := 2 


m := 2 


E : = l 
0 


q := 0 . .50 
v := 0 . .50 


Next calculate the 
profile:- 


relevant x and y Hermite polynomials for the chosen mode 


H 

x 


H 

y 



Finally calculate the intensity mode profile:- 


I(x,y) 


w 

r 2 2 t 

-» 

0 

- Lx + y J 


E ' H (x) • H (y) ■ exp 



0 W X y 

2 


- 

w 



The following variables 
chosen spot size, whilst 
routines. 


a M e i2 S S°S™ te<i With scal ^ n 9 the grid size to cover the 
is a dummy array reguired by the surface plotting 


x 

q 


q- 4 ■ w 

- 2 • w y 
50 


V- 4 • w 


50 


2 ■ w M 

q,V 
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The node profile 
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9 Appendix (3) 


Modifications to the LP-140 Pulsed COo Laser 

for Lidar Use 

Holger M. Jaenisch and Gary D. Spiers 
Center for Applied Optics 
University of Alabama in Huntsville 
Huntsville, Alabama 35899 


ABSTRACT 

The Pulse Systems Inc. Model LP-140 is a commercial 
pulsed CO 2 laser designed for marking and engraving. It 
is available with pulse energies in excess of 1 joule and 
repetition rates up to 7 Hz, thus making it a potential 
candidate for lidar applications. We document the char- 
acteristics of the LP-140 performance including power, 
temporal and spatial mode stability, chirp, and long ter- 
m operational characteristics. The laser can be made to 
function as a coherent lidar only if modified to improve 
its inherent characteristics. This paper addresses work in 
progress on the following modifications and their effect 
on performance: gas flow, optical resonator configuration 
and discharge supply modifications. 

INTRODUCTION 

The LP-140 is a transversely excited pulsed discharge 
device based on a sympathetic discharge 1 . The system 
is produced and marketed commercially for applications 
^uch as engraving. The LP-140 is designed and built to 
t he Pulse Systems Inc. (PSI) specifications 2 listed below. 

• DIMENSIONS: 0.35 x 0.35 x 1.2 m 3 (this dimension 
includes all optical components) 

• ENERGY PER PULSE: 3 J. (multi-mode with no 
frequency discrimination) 

• GAS CONSUMPTION: 0.3 Standard Cubic Feet Per 
Hour (0.14 1/min) 

• GAS MIX: 20% C0 2 ; 15% N 3 ; 65% He 

• INPUT POWER: 115V, 50-60 Hz. Q 5 amps 

• MAXIMUM PULSE RATE: 8 Hz ( Higher repetition 
rates available but at lower pulse energies) 

• PRESSURE: 40 torr 


• PULSE WIDTHS: 35 /is FWHM (100 ns to 200 /is 
available on special order) 

Our goal was to modify the LP-140 to meet the criteria 
outlined below. 

• Minimum overall size and compactness 

• Energy of I J per pulse 

• Maximum pulse length of 5 /is 

• Stable single mode operation 

• Interpulse frequency stability of < 1 MHz 

• Intrapulse frequency stability of < 0.2 MHz 

• Maximum pulse length of 5 /is 

• Maximum repetition rate of 10 Hz 

These requirements were derived from the requirements 
of the laser systems proposed for LAWS 3 (Laser Atmo- 
spheric Wind Sounder). The LP-140 was never consid- 
ered a contender for this program, however the relatively 
low cost, and simple, rugged design make it an attractive 
choice for small scale atmospheric studies if the laser per- 
formance can be improved. The LP-140 can also be used 
as a testbed for LAWS concepts. 

LASER DESCRIPTION 

The laser secured for this investigation differed from 
the standard commercial models on several key parame- 
ters. The active volume was extended from the «1.7 L of 
the standard design to «2.2 L. The standard gas mixture 
of C0 2 :N 2 :C0:He, 18:15:2:65 which produces long puls- 
es was modified to produce short pulses (< 10/is). The 
new mixture was designed by PSI using a Boltzmann dis- 
tribution code 4 and proprietary laser model 5 . The new 
mixture consists of C 0 2 :N 2 :C 0 :He, in the proportions 
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15.29:17 .33:1.99:65.39. The gas pressure is measured us- 
ing a Wallace & Tiernan model FA 160 gauge and the laser 
;s operated in the 30-50 torr pressure range. 

The standard LP-140 is configured for high power 
multi-mode operation. Our application requires mod- 
e discrimination, therefore the optical cavity design was 
changed to a hybrid, positive branch half-symmetric con- 
figuration using a Burleigh Instruments Model I- 1000 re- 
flection grating mounted at the Littrow angle to enable 
single longitudinal mode operation. The grating was 
also supplied with a Burleigh Instruments Model TS- 
100 translation stage, Model PZL-M Micrometer moun- 
t. Model PZL-015-00 15 /im PZL pusher, and a Model 
PZ-150-1 Amp/Driver-150 for feedback control and cavi- 
ty tuning. 

As presented in reference 2, each of the two gain arms 
has an effective volume of 4 x 4 x 70 cm 3 , comprising of 
tour equal discharge volumes, with an unsaturated gain 
t 0.031 cm’ 1 and a saturation intensity of 0.017 J cm’ 2 . 

The standard gas flow plumbing was engineered by PSI 
:o enhance operational lifetime. The commercial housing 
was maintained, but modified to allow the addition of 
the grating externally to the housing. The mechanical 
resonator structure was of a cast aluminum single piece 
construction. 

OPTICAL DESCRIPTION 


Our LP-140 folded optical cavity consists of a 6.35 cm 
diameter x 0.63 cm edge thickness meniscus ZnSe output 
coupler with a 5m radius of curvature on both surfaces 
.Figure 1). The A.R. coated convex side faces the dis- 
charge, in the center of which a square gold mirror had 
been flashed on with an area of 4 cm 2 . Two rectangular 
^old flat mirrors fold the optical path into the second ar- 
m of the optical cavity. Each is supported by a simple 
t hree-point mount. Finally the laser-tube cavity is sealed 
with a 10 m focal length, 6.35 cm diameter x 0.63 cm edge 
thickness A.R./A.R. plano-convex lens, with the convex 
side facing the discharge. The optical cavity is complet- 
ed by the grating placed « 10 cm beyond the final lens 
element. The optical cavity length is 238 cm. Since the 
4x4 cm 2 discharge region is square all optical elements 
are underfilled. The actual beam dimensions are (center 
obscured) 4x4 cm 2 . 

INITIAL CHARACTERIZATION 
Pulse Energy 

The initial pulse energy for the system was measured 
as approximately 1 joule per pulse at 1 Hz. The pulse en- 
ergy stability was poor and varied not only from shot to 



Figure 1: Optical layout of LP-140. 



PuIm ruftiKx) 

Figure 2; Pulse energy dependence upon pulse rate and 
laser temperature. 

shot but decreased to below a joule after only a few min- 
utes of operation. As the repetition rate was increased 
the average pulse power decreased as seen in Fig. 2. Dur- 
ing the test it was found that a minimum of 40 minutes 
was required for the operational system to reach thermal 
equilibrium and maximum stability. Also shown on the 
same figure is a similar result taken after the laser had 
reached thermal equilibrium. It can be seen that the out- 
put energy is consistently lower than immediately after 
startup. 

In the cool system the power varied by 10% from shot 
to shot, after the laser had warmed up this dropped to 
5%. 

Temporal Profile 

The temporal profile was measured using a PSI fast 
response pyroelectric detector, Model # UF-1 with a re- 
sponsivity of 10“ 7 V/W. To prevent detector damage the 
beam was attenuated before being focussed onto the de- 
tector. 

Figure 3 is an example of the initial temporal profile 
observed during the baseline characterization. It can be 
seen that there is significant mode beating during the 5 
/is pulse duration. 
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Figure 3: Pulse temporal profile prior to modification. 



Figure 4: Pulse spatial profile prior to modification. 


Spatial Profile 

The spatial profile was measured using a Spiricon beam 
profiler and the laser transverse mode pattern was found 
r o be highly unstable with a large degree of beam point- 
ing instability. Local hot spots in the beam distribution 
would vary on a shot to shot basis indicating possible 
transverse mode hopping. Figure 4 shows two represen- 
tative records. 

Frequency Chirp 

The intrapulse frequency chirp was measured by feed- 
ing the laser output into the existing Marshall Lidar 6 . 
The data was collected on the 10P20 line. Due to trans- 
verse mode instabilities indicated above the frequency 
chirp was considerable ( MHz) and varied inconsis- 
tently from shot to shot. 

SYSTEM MODIFICATIONS 
Gas Flow 

The degradation in output energy with both increased 



Figure 5: Average beam diameter dependence upon grat- 
ing extension. 


pulse rate and number of pulses tends to indicate insuffi- 
cient gas flow through the laser. This results in pressure, 
thermal, and possible dissociation mechanisms in the gas 
leading to the reduced output. During operation, the nee- 
dle on the pressure gauge momentarily rises by as much 
as 10 torr. This is indicative of the pressure wave which 
exists inside the cavity during the laser discharge. 

Inspection of the laser head showed that the gas in 
and out flow were limited by the orifice size of the flow 
tubing. The standard LP-140 laser flows gas through the 
discharge at a trickle rate of 0.5 L/min. The gas flow 
was modified by increasing the diameter of the input and 
output lines and reversing the direction of flow in the laser 
head. This enabled the flow rate to be increased to 2.5 
L/min. 

Cavity Length 

The cavity was lengthened by increasing the spacing 
between the grating and the lens, LI (Fig. 1). This 
enabled a more stable operation point to be obtained. 
Additionally the increased path length reduces spurious 
off-axis modes. The cavity was extended in 2.54 cm incre- 
ments. Previously, the grating would have to be adjusted 
every two to three shots to maintain a uniform laser out- 
put. For an extension of 53±5 cm the distribution became 
very uniform and stable. The optimum beam stability 
(Fig. 5) and mode pattern points did not coincide. The 
energy increased to approximately 1.25 joules (Fig. 6) for 
an extension of 51 cm but there was very slight mode in- 
stability characterized by the occasional appearance and 
disappearance of a mode hot-spot. 

At 54 cm extension the beam pointing instability was 
minimized. At 56 cm the energy decreased to 1.15 joules 
and the mode distribution was at its most stable point. 
The output quickly became chaotic as the cavity exten- 
sion decreased below 45 cm. There was also an accompa- 
nying energy drop to 0.75 joules for extensions below 25 
cm. As the cavity spacing was further decreased, the en- 
ergy once again increased. For extensions above 58 cm the 
power dropped off linearly. Mode and pointing instability 
never set in for cavity extensions out to 90 cm. Although 
stable, the output energy distribution collapsed to fill on- 
ly a single corner of the entire possible output pattern. 
This corner remained very stable in distribution but fluc- 
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Figure 6: Pulse energy dependence upon grating exten- 
sion. 

tuated in power. For extensions beyond 75 cm it became 
impossible to fill the output pattern despite adjustment 
to the grating. 

During these adjustments we found that rotating the 
grating by 90 0 decreased the output energy by 50-75%. 
We also rotated the output coupler by 45° and found that 
the spatial output pattern degraded severely. We feel this 
is due to diffraction and polarization effects caused by the 
sharp edges of the square reflector. We will address both 
of these phenomena in more detail in future work. 

When the lens LI was replaced with a flat uncoated 
NaCI window the laser operated with a stable output en- 
ergy and spatial distribution. The energy however de- 
creased by 50%. We feel this is accounted for by the lack 
of A.R. coatings and the marginal quality of the sub- 
strate available at that time. The increased diffraction 
losses resulting from the use of a flat window may prove 
ro be sufficient to cause stable operation without grating 
extension. 
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Figure 7: Pulse temporal profile alter modification. 

CHARACTERISTICS OF THE 
MODIFIED LASER 

Pulse Energy 

The pulse energy was immediately observed to have 
increased slightly and stabilized. At the most stable point 
for pulse to pulse repeatability, the energy is equal to the 
original performance of 1 J /pulse. When the cavity length 
is slightly modified from this point by an extension of 2-5 
cm, the energy increases by 10—15%. The spatial profile 
however is not quite as stable. 

The optimum point of operation Is better described as 
a window about 11 cm wide peaked at a cavity length 
extended by 56 cm. 

Temporal Profile 

The temporal profile at the peak cavity length improves 
substantially (Fig. 7). The pulse length not only short- 
ens, but the tendency for multi-pulsing is also reduced. 


Matrix Method LP-140 Model 

A MathCad 1 ^ version 2.5 document is presented in 
Appendix A and may be used without alteration. The 
matrix technique used is taken from Siegman 10 . 

From this document, it can be seen that the extended 
cavity has a magnification of 2.0 and an equivalent Fresnel 
number of 1.7. On extending the cavity by 0.56 m and by 
changing the value of D2 in the document, the magnifica- 
tion is increased to 2.1 and the equivalent Fresnel number 
changes to 1.6. This implies that the improved mode sta- 
bility is due to combination of increased diffraction losses 
and slightly improved Fresnel number. 


At the original cavity length, inconsistent pulse lengths 
up to 12 fss were not uncommon with triple spikes in 
the pulse. By extending the cavity, these pulses have 
been reduced to 5-8 /is with a single weil defined spike. 
It is important to note that at all cavity spacings, the 
gas mixture and pressure also play an important role in 
defining the temporal profile. 

Spatial Profile 

The spatial profile at the optimal cavity spacing demon- 
strates the most significant improvement (Fig. 8). Out- 
side of the 11 cm window the profile is chaotic and exhibit- 
s strong beam pointing instabilities. Within the window 
the temporal profile fills out and the beam settles down 
to a high degree of stability. 

It is difficult to quantify the exact nature of the spa^ 
tiai mode. The center obscuration creates a square out- 
put beam with & square hole in the center. The f&r-field 


54 









Figure 8: Pulse spatial profile after modification. 
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diffraction pattern is non-uniform and may not be ade- 
quately described as gaussian or top-hat. Both the near Figure 9: Frequency chirp, 

and far-field patterns display large amounts of diffraction 
ripple. 

Frequency Chirp 


The chirp data was taken for pressures of 30-50 torr 
and was found to be stable and repeatable with a de- 
pendence on pressure. The higher pressure range from 
40-50 torr provided the best results. The frequency chirp 
is shown in Figs. 9 and 10. The initial large downsweep 
in frequency can be attributed to the decaying electron 
density in the discharge' while the upswing is the conven- 
tional t 3 time dependence 8 . It can be seen that the pulse 
is approximately 13 /is long compared to the apparent 5 
lis in the pulse temporal profile. This apparent difference 
is caused by the heterodyne technique which extracts the Frequency 
long, very low energy tail of the pulse from noise. Over c ^ r P l kHz) 
t his total pulse the frequency chirp is approximately 3 
MHz. In practice the latter portion of the pulse is in- 
significant to the application and can be ignored. It can 
be seen from Fig. 10 that there is a 5 /is long portion of 
the pulse with approximately 250 kHz of chirp which is 
acceptable for our application. It should be possible to 
eliminate the leading and tail edge of the pulse to provide 
a slice which would contain 70% of the pulse energy. 



Time (|is) 


Discharge Characteristics 


Each of the discharges is driven by a separate cir- 
cuit (Fig. 1 1), although they are all switched by a com- 
mon EG LG HY1102 thyratron. The discharge voltage 
pulse was measured using a Tektronics P6015 high volt- 
age probe and the current pulse using a Pearson Elec- 
tronics Model 110 induction coil. Both were connected 
to a LeCroy 9450 digitizing oscilloscope which download- 
ed the signal pulses to an IBM PC compatible computer 
for analysis. Typical current and voltage pulses for one 


Figure 10: Frequency chirp. 
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Figure 11: Sympathetic discharge circuit. 
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Figure 12: Voltage pulse. 


of the electrodes are shown in Figs. 12 and 13. It can 
be seen that there is considerable ringing on the voltage 
pulse, which is due to a mismatch in the impedance of 
the discharge circuit and the laser discharge. 

The electrical energy into the discharge can be obtained 
by integrating the voltage and current pulses (Fig. 14). 
This shows that the energy deposited into a single dis- 
charge is %3.5 J. However, only the energy deposited into 
the CO2(001 ) and the 1^3(1/ =1—8) vibrational energy 
levels is useful for lasing action. Lowke et al. 9 have shown 
that the optimum excitation of these vibrational levels oc- 
curs for a reduced electric field or E/N of « 1 — 3 • 10” 16 
Vcm 2 . The value of E/N was derived from the discharge 
voltage pulse form for each digitization interval together 
with the energy input in that interval. This enabled the 
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Figure 13: Current pulse. 



Figure 14: Energy over time. 
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Figure 15: Energy input. 


input energy into the discharge to be evaluated as a func- 
tion of E/N for E/N intervals of 1 10” 16 Vcm 2 (Fig. 15). 
It can be seen that most of the energy is deposited in the 
4 — 6 -10“ 16 Vcm 2 . This results in the low laser efficiency 
of =s4% compared to as 12% for the laser in its original 
multimode long pulse configuration. This means that the 
discharge circuit requires optimizing for the new gas mix- 
ture. The evolution of the discharge impedance with time 
was found by dividing the voltage by the current (Fig. 
16). It can be seen that the discharge impedance during 
the bulk of the energy input vanes between 5 ft and 20 
ft. This enables the discharge circuit to be optimized. 



10 10 413 


Figure 16: Discharge impedance. 


56 



